Massless metric preheating 
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Can super-Hubble metric perturbations be amplified exponentially during preheating ? Yes. An 
analytical existence proof is provided by exploiting the conformal properties of massless inflationary 
models. The traditional conserved quantity £ is non-conserved in many regions of parameter space. 
We include backreaction through the homogeneous parts of the inflaton and preheating fields and 
discuss the role of initial conditions on the post-preheating power-spectrum. Maximum field vari- 
ances are strongly underestimated if metric perturbations are ignored. We illustrate this in the case 
of strong self-interaction of the decay products. Without metric perturbations, preheating in this 
case is very inefficient. However, metric perturbations increase the maximum field variances and 
give alternative channels for the resonance to proceed. This implies that metric perturbations can 
have a large impact on calculations of relic abundances of particles produced during preheating. 
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I. INTRODUCTION 

Reheating is a crucial element of any inflationary cos- 
mology. It endows the nascent universe with the high 
temperatures and large entropy needed for its subsequent 
evolution. Reheating changes the equation of state from 
the near de Sitter p ~ — p, to the radiation-dominated 
form p = This leads to a large amplification of the 
power spectrum of super-Hubble metric perturbations re- 
sulting in a final value Pk oc (f>~ 1 \k=aH , where <f> is the 
slowly-rolling inflaton field [0J|]. This we call the 'old' 
theory of cosmological reheating. It is a theory where 
coherence of the inflaton, and the fact that the Hubble 
scale, at the end of inflation is dwarfed by the vastly 

larger particle horizon dn, play no role. 

In contrast, the super-Hubble coherence || of the infla- 
ton condensate forms the foundation of preheating Q . It 
is an intrinsically non-equilibrium era || in which stim- 
ulated, coherent processes can totally dominate single- 
body, perturbative decays. The resulting phase of vio- 
lent particle production leaves as its legacy a non-thermal 
spectrum of particles. The subsequent backreaction and 
rescattering [g| of these fields leads to turbulence and 
thermalisation 0. 

While many different models have been proposed, the 
basic mechanism of parametric resonant amplification of 
quantum fluctuations is very robust and bears many sim- 
ilarities to standard quantum field theory in strong ex- 
ternal fields H§|. Given this robustness, a natural ques- 
tion is whether these time-dependent mass effects have 
implications for the standard predictions of inflationary 
cosmology 0. Primary among these predictions is the 



spectrum of metric perturbations, which will provide pre- 
cision tests of the inflationary paradigm. It is therefore 
of great interest to understand these issues in detail. 

That parametric resonance due to oscillating scalar 
fields might have an interesting impact on metric per- 
turbations was first noted by Kodama and Sasaki [||. 
The first works addressing this issue within the context 
of preheating however claimed that there would be little 
effect on metric perturbation evolution on super-Hubble 
scales, and in particular that C (see eq. fl0| ), would remain 
conserved |^0| , pT[ . These conclusions stem from consid- 
ering small resonance parameters q < 0(1) such as oc- 
curs in the single field \<fi A model jlj . While these models 
can yield a rich spectrum of sub-Hubble physics |t2| they 
cannot be considered generic models of preheating, which 
usually requires q 3> 1 in order to be efficient Q . 

Detailed simulations |l4|] have presented strong evi- 
dence in support of the intuitive arguments Q for the 
existence of robust resonances on super-Hubble (a.k.a. 
"super-horizon" ) scales |Q . Such amplification does not 
violate causality |l3|,|l4|. Nevertheless, an analytical ex- 
istence proof of the possibility of exponential growth of 
super-Hubble metric perturbations remained outstand- 
ing. This Letter aims to provide this missing link. We 
will show that preheating can in principle cause signifi- 
cant departures from the old theory of metric perturba- 
tion evolution through reheating. 

The failure of the old theory is most glaringly seen 



*It is important to distinguish between metric amplifica- 
tion due (1) purely to resonances which stem from reversible 



changes to effective masses and (2) to thermalisation which 
involves irreversible change in the average equation of state 
and large entropy production. 

^Such as the fact that gravity has negative specific heat: 
inhomogeneity is the natural end-point of gravitational 
evolution. 
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in the violent, essentially exponential^], growth of the 
Bardeen potential (. In the old theory, £ is exactly 
conserved in the adiabatic, long-wavelength limit. This 
allowed one to match the scalar metric perturbation $ 
across the desert separating the Grand Unified scale and 
photon decoupling where temperature anisotropies in the 
CMB were formed on large angles with AT/T ~ i$ 
(fi = l,A = 0). 

In section (|l|) we set out the general class of models we 
consider together with the general, conformally rescaled 
multi-field equations. Section (III) generalises the work 
of |ij]|l8| to include metric perturbations. This provides 
the basic analytical results of our paper. Section (IV) 



discusses the appropriate initial conditions for the x ne ld 
at the start of preheating while section ((v|) discusses the 
X power spectrum at the start of preheating, one of the 
key factors in evaluating the impact of preheating on the 
CMB. Finally section ( |v| ) discusses the effects of back- 
reaction and self-interaction among the decay products. 
We use natural units where G = 1. 



II. THE GENERAL CLASS 

Our general potential for the class of massless scalar 
fields (fit, interacting with couplings gijhU will be: 



JV 

This includes many interesting sub-classes such as the 



archetypal X4> 4 +g 2 4' 2 X 2 model considered in section (III) 



We choose this form since it is the most general one con- 
sistent with conformal invariance [j]. The fluctuations in 
these fields generate scalar metric perturbations |**| which, 
in the longitudinal gauge, are encoded in the perturbed 
metric 



ds 2 = a(i]) 2 [(1 - 2$) di] 2 - (1 + 2*) CbfdXi] 



(2) 



where dr\ = dt/a is the conformal time, a{rf) is the scale 
factor, and $ = if? since the anisotropic stress of the 
system vanishes to linear order |l4|] . It proves very useful 
to re-scale all the fields by a(rj). We adopt the convention 
that F = aF for any field F. The spatially homogeneous 
parts of the fields satisfy the background Klein-Gordon 
equations: 



(p'l - ^ + Vfc = , (3) 

where / = d/dr/, Vm i = dV/d(pi and V((p) = a 4 V((p) is 
the conformal potential, given by 

V(<fi) = £i,j,i,m 9ijlm <Pi t Pj<Pl<Pm ■ (4) 

The perturbed, multi-field, Einstein equations are: 
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(5) 
(6) 



We now reformulate this system in terms of the rescaled 
Sasaki-Mukhanov variables: 



a' a 



(7) 



which satisfy the equations: 
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The constraint equation (0) can be recast as: 



k 



2H- a - 



4-7T 

<z> k = -J2(^-n?)Qi, (8) 



where Ti = a' /a. Note that when a" — 0, the LHS of 
these equations are expansion-invariant (EI) while the 
non-EI coefficients on the RHS are: 



M lJ = 
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Finally, the traditional "conserved" quantity £, rescaled 
by a, is given by: 



(k = 



a 



(10) 



Equations (|7|-|Tc|) will provide the foundation for our an- 
alytical and numerical results. 



"•"By essentially exponential we mean exponential modulo 
sub-dominant power-law factors; i.e. e M /t a , a < 1. 

§We neglect the conformal coupling to the curvature since 
the associated resonance is weak for £ = 1/6 

**We neglect vector and tensor perturbations |^(| and con- 
sider a flat FLRW background. See plpl for aspects of the 
non-flat case. 



III. AN ARCHETYPAL EXAMPLE 

Within the general class of potentials (Q), a key, well- 
studied example is the massless \<fi model coupled quar- 
tically to a massless \ jl^,[l8]]. The conformal potential 
in this case is 
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(ii) 



We switch to dimensionless time x = a/A^o 7 ? where 
ipo — 0.3 is the value of 4> at the start of preheating 
corresponding to xo — 2.44 In this model, the time- 
averaged equation of state is that of radiation and the 
scale factor obeys: 



(2^)1/2 ^=(f)V2 ^ 



This implies that a" 
satisfies: 



0. Define (p 



v 3 + ik^ = °- 



(12) 

4>/y>o- Then if 
(13) 



which during the first stage of preheating when x — 0, 
has the oscillating solution ||[l7| >f{x) = cn(x—x , 1/V%) 
where an is the Jacobi elliptic cosine, while 



X 
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0. 



The Ql satisfy: 



(14) 



weak, super- Hubble, growing mode to Q\. The system is 
decoupled, exact Floquet theory can be applied to the Q\ 
equation and we are guaranteed solutions growing expo- 
nentially as e^ K ' x for values of n and g 2 /A inside resonance 
bands; see Fig. ([[]). Upon neglect of M vv , Floquet the- 
ory can also be applied to the Qf equation. 
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FIG. 1. Floquet index, /j, k , instability chart. White repre- 
0, darker greys represent increasing fi K (left). fj> K 



sents fj, K 

vs k 2 = k 2 /\ip 2 for g 2 /\ = 2 (right). Both plots are for the 
generalized Lame equation given by Eq. (|l(j) with the RHS 
set to zero and <f> = cn(x, l/\/2). 
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(15) 



(16) 



where k = k 2 j\ip^. These equations split naturally 
into terms which are expansion- invariant (EI) , and those 
which are non-EI which come from the M u in eq. (j|) 
and which depend explicitly on a. We now consider these 
equations with and without x backreaction. 



1. Case 1: Reduction to generalized Lame form 

Let us first study the equations (p^lq) without any 
backreaction i.e. using (p = cn(a;,l/v2). This elimi- 
nates all terms on the RHS except that containing M w . 
We verified numerically that this term provides only a 



^By backreaction we mean including terms involving prod- 
ucts of the background \ field and a genuine, inhomogeneous, 
first order variable such as Q\. While \ ls homogeneous, 
so that the product is ostensibly first order, % i s vanishingly 
small during inflation. Hence it is numerically irrelevant at 
the start of preheating. Once x grows large it cannot be ne- 
glected and leads to the termination of the resonance (see Fig. 
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2. The gauge-invariant metric perturbations $ and 
£ for the super-Hubble mode k 2 — 10 -20 and p 2 /A = 2. 



Note that $fc grows with Floquet index /i ~ 0.357, much 
larger than the maximum possible in preheating neglecting 
metric perturbations which is bounded for all values of g 2 /X 
to be less than or equal to 0.238. Inset: for g 2 /\ = 2 
exponentially growing) and g 2 = (£k constant). 

In fact the equations are of generalized Lame form in 
Minkowski space, as studied in depth by Greene et al 
JItJ and Kaiser 18 1. The Floquet index /i K is shown in 
Fig. (0) as a function of g 2 /\ and k 2 . For g 2 /X = 2 
it reaches its maximum value of /i K ~ 0.238 at k 2 = 0, 
which coincides with the global maximum of the Floquet 
index. Further, there exist an infinite number of reso- 
nance bands in g 2 /X for which /i K =o > 0, showing that for 
large regions of parameter space there is essentially ex- 
ponential growth of super-Hubble metric perturbations. 
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This provides the first conclusive analytical proof of this 
fact, confirming previous numerical work |r^ , |l4| ]. 

In Fig. §) we plot £ fe , Ql and Q\ for k 2 = 10" 20 and 
g 2 /X = 2. Note the perfect exponential growth of Q\, 
while Qf, simply oscillates with constant amplitude since 
it corresponds to g 2 /X = 3 for which there is only a sub- 
Hubble resonance band. The Floquet index for is as 
expected from the Lame analysis, while the Floquet index 
for <i>fe is much larger, around 0.357. This is expected 
from the constraint equation (||) which shows that is 
sensitive to the product xQk j with both \ and Q\ growing 
exponentially. 

Expressed in dimensionless time and neglecting the 
a" ~ term in eq. (|To|), f becomes 



(17) 



which shows that C,k also grows exponentially, sharing the 
same Floquet index as before backreaction shuts off 
the resonance. This is evident in the inset of Fig. (||) 
where we show £ for g 2 /X = (no growth) and g 2 /X = 2 
(exponential growth). 

Finally, note that the effective resonance parameter 
q e ff — g 2 /X can be as small as q e ff — 2 for the Floquet 
index to reach its global maximum and for super-Hubble 
preheating to be strong. This means that the x is 
not necessarily suppressed during inflation p3[ . This is 
discussed in detail in section (IV). The suppression also 



does not apply in wide classes of other preheating models 



2. Case 2: Backreaction included 

We now consider the full equations ([l3] - |l6|) without 
approximation. Fig. ([|) shows the results of simulations 
including backreaction on the inflaton. For x < 50, the 
analysis in case 1, presented above, applies. However, 
when lnx ~ 0(1), the backreaction in Eq. ( p^| ) becomes 
important, changing the amplitude and frequency of (p 
oscillations as shown in the inset of Fig. (H) . This pushes 
X into a region of narrow resonance with ~ 0.05. 

Importantly, however, the metric perturbations do not 
stop growing once x moves into the region of narrow res- 
onance. The couplings between Eqs. (j^) and ( |l6| ) cause 
Q\ to begin growing exponentially and synchronisation 
occurs: for x > 50 the perturbations and Q\ all 

grow with the same, large, Floquet index. This shows 
that any analysis based on the decoupled equations one 
finds in the absence of metric perturbations must be con- 
sidered with caution, since it would (1) underestimate the 
Floquet indices and (2) incorrectly predict no growth of 
Qfc and 5tpk for k 2 « 1. From Eq. (R) one sees that the 
field fluctuations 6<pk grow with the same, large Floquet 
index as the metric perturbations. 



Backreaction finally ends the growth of the metric per- 
turbations around x ~ 135 when the amplitude of <j) 
drops and the inflaton effective mass dramatically in- 
creases (see Fig. 0). 
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FIG. 3. The field \ an d the gauge-invariant metric pertur 
bations Q\ and Q% for k 2 = 10 
the onset of backreaction in x at i 
sation in the metric perturbations for x > 50. Inset: The 
evolution of the inflaton condensate <f>: note the drop in am- 
plitude and increase in oscillation frequency at x ~ 50 and 
135. 



IV. PREHEATING INITIAL CONDITIONS 

The key point that arises from having resonant growth 
of the super-Hubble metric modes with k <C 1 for 1 < 
g 2 /X < 3 is its relation to initial conditions at the start 
of preheating. Convolving these two will yield the post- 
preheating power spectrum. 

The x effective mass is given by 



°x,eff 



2i2 

9 <i> 



(18) 



Ignoring the metric perturbations during inflation where 
they are small, one can solve the Klein-Gordon equation 
in de Sitter spacetime for the x modes in terms of Hankel 
functions |E3l: 



Xk - c ly ftHi 2 \kr,) + C2i/vH^ ) (kr ] ) , 



(19) 



where c\ — C2 = corresponds to the adiabatic 

vacuum state and the order, is, of the Hankel function is 
the root of 



H 2 ^ 



(20) 



Here £ is the non-minimal coupling of the x field to the 
curvature, which we have put to zero throughout so far. 
Asymptotically, as kr\ — > 0, the Hankel functions are pa] 



H 



(1,2) 



(M -> -r(i/) 



(21) 
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so that if Re(V) > 0, the mode functions are infra-red 
divergent. However, when Re(z/) = 0, it is simple to 
show that the mode functions are finite as kr\ — > and 
are essentially /c-independent. Further, if Re(z/) = 0, 
the modes decay away as a~ 3 / 2 , that is, exponentially in 
proper time. This is the origin of the claim that preheat- 
ing has no effect on super-Hubble metric perturbations 

Now, in our case, since x — minimises the potential 
at fixed <fi, the x field has no potential energy during 
inflation and H 2 ~ 27rA0 4 /3. Hence when £ = 



9 
4 



3 5 2 



27TA0 2 



(22) 



hence for g 2 /X = 2, v 2 < when < 2/V3tt ~ 0.7M pl . 
Now since preheating starts at around 4> ~ 0.3M p i, v 2 > 
during most of inflation. This case differs radically from 
the case g 2 /X S> 1, where v is complex many e-foldings 
from the end of inflation. Hence, while the suppression 
of x modes is very strong when g 2 /A ^> 1 p^| , for g 2 /X 6 
[1, 3) there is little suppresion of cosmological x modes. 

Before discussing the post-preheating power spectrum 
we make two points. Firstly, our potential, Eq. (|i"T|), is 
not meant to be generic, and hence specific results should 
be taken as proofs of principle rather than as generic. In 
particular, the suppression of x modes when g 2 /X ^ 1 is 
absent in preheating with more general potentials 

Secondly, the suppression that occurs for g 2 / \ 2> 1 and 
X = (minimally coupling), does not necessarily persist 
when £ 7^ 0. From Eq. (|5|) one sees that v 2 > is 
compatible with to 2 ^> H 2 (and hence little or no sup- 
pression) for large negative couplings, £, to the curvature. 

In fact, for sufficiently large negative £, v > 3/2 and 
the resulting power spectrum has most power at large 
scales - a red tilt - and field variances receive their dom- 
inant contribution from the super-Hubble modes 
Large non-minimal couplings in fact lead to their own 
strong resonances due to the oscillation of the Ricci scalar 
during preheating p6 26 1. 



V. THE x POWER SPECTRUM 

While our main aim in this paper has been to pro- 
vide analytical evidence for the existence of super-Hubble 
resonances, we turn here to the question of the post- 
preheating power spectrum. 

The power spectrum for an arbitrary field X is defined 

by 



(23) 



where (•) denotes ensemble averaging (equivalently vol- 
ume averaging if X is ergodic). From Eqs. ([L9,Elf), the 



power-spectrum for x a t the end of inflation in the limit 
krj ^ Q is thus given by: 



P X (M) 



= &" r ^ (t 



-2v 



(24) 



This clarifies the situation immediately. For v complex, 
the power spectrum is steep oc fc 3 . However, for v real, 
the power spectrum is closer to Harrison-Zel'dovich with 
spectral index 



3-2z/, 



(25) 



which tends to zero as v — ► 3/2 (<fi — ► oo). Now it is true 
that v becomes complex at the end of inflation, but since 
we are interested in the k — ► part of the spectrum, the 
appropriate value of v is the one at early stages: <f> > 
2M p i, for which v > v2j leaving only a mild blue tilt. 

The power spectrum ( p4| ) is then modified by preheat- 
ing in a rather trivial way: it is multiplied by a factor 
~ e 2 ^ Ax , where Ax is the time spent in resonance. For 
k 2 <C 1, modes, this transfer function is almost indepen- 
dent of k since the Floquet index only changes signifi- 
cantly over Ak ~ 0(1). 

The above discussion paves the way for direct evalu- 
ation of the impact of preheating on large scales. The 
obvious question remains: "how much of the g 2 /X phase 
space will leave an imprint on the CMB?" In this model 
we may investigate this question analytically without re- 
lying on numerical studies. 



A. Dividing up the g 2 /\ parameter space 

For cosmological implications we are typically only in- 
terested in the k 2 < 10 -50 section of the n 2 — g 2 /X chart, 
Fig. (0). It is relatively easy to show that the nodes of 
the resonance bands (where (J. K =o becomes non-zero) are 
given by g 2 /X = n(n + l)/2, with n natural. It is then 
easy to show that asymptotically for large g 2 /A, half of 
all values of g 2 / X lie in resonance bands, and half lie in 
stable regions. 

However, if £ = 0, only for a much smaller measure 
of g 2 /X values do the cosmological modes dominate the 
backreaction integral, immediately ruling out those val- 
ues of g 2 IX. This clinical division of the g 2 jX parameter 
space is due to the conformal invariance of our model. 
For generic models with massive fields the modes will 
move through the maxima of the resonance bands and 
the issue of which couplings are ruled out is more subtle, 
and probably needs to be done numerically. 

On the other hand, modes near g 2 /X = 1 or 3, will 
experience super-Hubble growth and are also little sup- 
pressed during inflation. However they are unlikely 
to go nonlinear before backreaction becomes important 
(since sub-Hubble modes grow the quickest - see Fig [l]) . 
Whether these regions of the resonance band are ruled 
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out is therefore more subtle, but they are expected to 
leave an imprint on the CMB. 

For g 2 1 A 3> 1 (£ = 0), the \ spectrum is suppressed 
at small k and although super-Hubble modes experience 
resonant growth, backreaction may well shut off the reso- 
nance before they become important and affect the power 
spectrum in any significant way. Equivalently n ~ 3 in 
Eq. @ |§. 

The g 2 /X parameter space is thus divided into four 
sub-regions: (1) regions where backreaction comes from 
cosmological modes going nonlinear. This is a small 
set in the conformal model and can be ruled out di- 
rectly. They form narrow regions around the maxima 
of the Floquet index at g 2 /A = 2, 8, .... (II) regions with 
< fx K - < ^maximum-, i.e. in the wings of the res- 
onance bands. The Q x variance is dominated by sub- 
Hubble modes but they can nevertheless leave an im- 
print on the CMB if g 2 /X < 10. (Ill) Regions with 
Hk=o > but g 2 /X > 10 where suppression during infla- 
tion is strong and little effect is seen. (IV) Stable bands 
where /ifc=o = and backreaction is irrelevant. 

VI. x SELF-INTERACTION 

We now consider the dual issue of the possible effects 
of the rapidly growing metric perturbations on particle 
production. 

A result from early studies of preheating, believed to be 
rather robust, was that \ self- interaction lead to a rapid 
shut-off of the resonance and inefficient decay of the in- 
flaton |22j . This is easy to understand in the absence of 
metric perturbations since the self-interaction acts like an 
effective mass m 2 , e ^ ~ 3A x (x 2 ), which pushes all modes 
out of the dominant first resonance band and essentially 
ends the resonance. This is evident in the inset of Fig. (Q) 
which shows $ in a model with A x x 4 /4 self-interaction 
added to the potential (|ll| ), with A x = 10~ 8 A for case 
1, i.e. the equations ( p^|Jl6^ with all terms on the RHS's 
set to zero and (p = cn(x, l/y/2). Even with this tiny 
self-interaction, backreaction appears to be extremely ef- 
ficient at ending the resonant growth of super-Hubble 
fluctuations. 

This, however, is an artifact of neglecting metric per- 
turbations which manifests itself as the couplings be- 
tween the equations (|l5| , |l6|) . The main graphs in Fig. 
(||) show a very different behaviour, even when a much 
larger self-interaction, A x = 10 2 A, is used to emphasise 
the discrepancies with the old picture. The metric and 
field fluctuations continue to grow exponentially long af- 
ter the x fi^d has ended its growth. How is this to be 
understood ? 

The growth of x forces it out of the resonance band 
and causes it to stop growing before it has a chance to 
affect the inflaton evolution through Eq. (|i"3|). This is 
in complete contrast to the case in Fig. (0). Since it is 



the amplitude and frequency of inflaton oscillations that 
predominantly controls the resonant growth of field and 
metric fluctuations, the large x self-interaction actually 
aids the growth of metric perturbations in the sense that 
it removes the means of backreaction on the inflaton con- 
densate 0. 
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FIG. 4. Evolution of metric perturbations with \ 
self-interaction: (a) $ fc (b) Q* (c) Q£ and the x field, for 
A x = 10 2 A, k 2 = 1CT 20 and g 2 /X = 2. While self-interaction 
stops x from growing, it fails to stop exponential growth of 
the metric perturbations. Inset: <3>fc evolution with all terms 
on the RHS of Eq's ([L5[[l|) set to zero, and A x = 10" 8 A. The 
apparent shut-off of the resonance is an artifact due to neglect 
of backreaction and coupling between and Q?. 

This provides an elegant example of both the qualita- 
tive and quantitative importance of including metric per- 
turbations in preheating. This is true not only when one 
is interested in understanding the spectrum of perturba- 
tions for the CMB, but also when calculating relic abun- 
dances of particles produced during preheating. Rele- 
vant examples are the abundances of gravitinos, GUT 
bosons and super-heavy dark matter whose decays may 
contribute to the observed ultra-high energy cosmic ray 
flux (27). 

VII. CONCLUSIONS 

In conclusion, we have attempted to convince the 
reader of a number of controversial possibilities: (1) 
gauge-invariant, super-Hubble scalar metric perturba- 
tions can grow exponentially, modulo power-law be- 
haviour, in an expanding universe. We did this by re- 
ducing the gauge- invariant equations to generalised Lame 
form. (2) The quantity £, often used to transfer the power 



**Using (x 2 ) instead of x 2 in Eq. (Q will cause the field 
fluctuations in Fig. (^) to stop growing at some large, but 
finite, value. However, our main point is to stress the new 
channels through which the resonance can proceed when met- 
ric perturbations are included. 
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spectrum of fluctuations to decoupling by virtue of its 
constancy, may be useless in general during preheating 
due to its exponential amplification. (3) Metric pertur- 
bations can significantly alter estimates for the maximum 
variances of field fluctuations needed for accurate predic- 
tion of relic particle abundances. 

This last point is particularly transparent when the in- 
flaton decay products exhibit strong self-interaction. In 
the absence of metric perturbations the self-interaction 
stops preheating almost completely When metric 

perturbations are included however, the resonances pro- 
ceed through new channels and can be strong even at 
large self-interaction, greatly enhancing final field vari- 
ances and number densities. These issues will be dis- 
cussed in depth in [fisf] . 

The issue of the post-preheating x power spectrum has 
been discussed as a function of the coupling g 2 /X. We 
show that while there are an infinite number of super- 
Hubble resonance bands, only for a small range of values, 
g 2 /A < 10, can one expect modifications of standard in- 
flationary predictions for cosmology since these values are 
partially or completely protected from strong suppression 
during inflation. In more general, non-conformal, models 
the situation will be more complex as modes are pulled 
through resonance bands by the expansion. Crucially the 
mode suppression mechanism may also be absent p5|. 

To end we discuss an issue little studied to date. Elim- 
inating all other perturbation variables, the $ equation 
of motion typically has singular coefficients oc [EH. 
Indeed, this was the main original motivation for using 
the Sasaki-Mukhanov variables |T^,|ll|. Insight into <f> 
evolution is provided by the duality between the Klein- 
Gordon and I-d spatial Schrodinger equations exploited 
in ]2l[ | . In this approach the resonance band structure in 
fc-space becomes dual to the spectrum of the correspond- 
ing Schrodinger equation. Now consider the dimension- 
less Scarf Hamiltonian: 



dz 2 



A 



■ 2 

sin z 



(26) 



whose potential is singular periodically when sinz = 0. 
This Hamiltonian is known to have a spectrum with band 
structure for —1/4 < A < and discrete eigenvalues 
otherwise S|. On the complement of the spectrum, 
R — cr(iJ), the corresponding Floquet index is real and 
positive pif| . This explains why the behaviour of $fc can 
be qualitatively similar to that of Q\ and Q^, which obey 
non-singular equations. 
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